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Abstract 

The asymptotic behavior of the optical potential, describing elastic scatter- 
ing of a charged particle a off a bound state of two charged, or one charged 



' and one neutral, particles at small momentum transfer A a or equivalently 

00 ' 

at large intercluster distance p a , is investigated within the framework of 

CO 

o 
>n 

Q\ • tion that occurs in the exact expression for the optical potential, a recently 



the exact three-body theory. For the three-charged-particle Green func- 



derived expression, which is appropriate for the asymptotic region under 
consideration, is used. We find that for arbitrary values of the energy pa- 
rameter the non-static part of the optical potential behaves for A a — ► as 
C\A a + o(A a ). From this we derive for the Fourier transform of its on- 

5— i 

shell restriction for p a — > oo the behavior — a/2p* + o (l/p„), i.e., dipole or 
quadrupole terms do not occur in the coordinate-space asymptotics. This 
result corroborates the standard one, which is obtained by perturbative 
methods. The general, energy-dependent expression for the dynamic po- 
larisability C\ is derived; on the energy shell it reduces to the conventional 
polarisability a which is independent of the energy. We emphasize that the 
present derivation is non-perturbative, i.e., it does not make use of adiabatic 
or similar approximations, and is valid for energies below as well as above 
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the three-body dissociation threshold. 
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I. INTRODUCTION 



The incorporation of the long-ranged Coulomb interaction in the theoretical descrip- 
tion of reactions between charged composite particles constitutes a problem of long- 
standing interest, due to its great importance in atomic and nuclear physics. As a proto- 
type we consider the elastic scattering of a charged particle off a cluster composed of two 
charged (or one charged and one neutral) particles. 

Formally, such a reaction can be described in terms of a single- channel Lippmann- 
Schwinger-type integral equation for the elastic scattering amplitude, in which the so- 
called optical potential occurs (see, e.g., Although the latter is structurally too 
complicated for practical purposes (besides being non-local, energy-dependent, and com- 
plex above the threshold for the opening of the lowest inelastic channel, it contains the 
three-body Green function which is actually the object whose calculation is attempted), 
its very existence is of considerable significance for many purposes. 

An alternative formulation is based on the exact three-particle equations ||, suit- 
ably generalized to accommodate the long-ranged Coulomb potentials either in coor- 
dinate f|], or in momentum space [0-0. The latter approach leads to coupled multi- 
channel Lippmann-Schwinger-type equations whose solution yield simultaneously all two- 
fragment, i.e. the (in-)elastic and rearrangement, amplitudes (and with an additional 
quadrature also the dissociation amplitudes). For three charged particles, the effective 
potentials occurring therein are again too complicated to be used presently for practical 
calculations (except with drastic approximations). But, if only two of the three parti- 
cles are charged, some version of these equations becomes manageable, and has in 
fact already been applied successfully to the calculation of elastic scattering and breakup 
observables in the proton-deuteron system |L0|L1 |. 



Due to the differences in the physical pictures on which these two formulations are 
based, the expressions for the optical potential and for the effective potential as used in 



conventional applications of the three-body theory, differ appreciably (indeed, in ref. JT2 
it has been shown that the former can be derived within the framework of the three- 
body theory only under very special and unusual assumptions). It is therefore of general 



3 



interest to gain insight into such properties of effective potentials, which do not depend on 
the specific details of the formulation, and are therefore common to both of them. One 
important example is the behavior for large separation of the colliding particles which 
should be similar in both cases, in contrast to the short-range behavior which in general 
will be rather different. 

In the present investigation we concentrate on the optical potential. It is known to 
simplify considerably when the distance between the colliding particles goes to infinity, 
or equivalently when the momentum transfer goes to zero. In fact, its large-distance 
behavior has been extracted long time ago by various perturbative methods below |13| , |T4|1 , 
and recently also above |15| (cf. also [0,0]) the dissociation threshold. The result is well- 
known: besides the static potential which comprises the multipole contributions arising 
from the charge distribution of the composite particle, the first non-vanishing term of 
longest range is the (local) polarisation potential (proportional to the inverse fourth power 
of the distance between the two colliding bodies). But in none of these investigations the 
reliability of the approximate, or the convergence of the perturbative, treatment has been 
studied. 

However, precise knowledge of the asymptotic behavior of the optical potential is of 
great interest, for many reasons: not only does it constrain the construction of (the long- 
range part of) model optical potentials for practical applications; but it also provides 
answers to crucial questions like, e.g., what kind of effective-range expansion be adequate 
or what rate of convergence of partial-wave sums is to be expected. 

For these reasons, there have recently appeared several attempts to provide non- 
perturbative derivations based on three-body theory, both in coordinate |18| and in mo- 
mentum [p!9| -|23[| space. Indeed, it was found that the above mentioned, approximately 
derived long-range behavior could be recovered. However, not only are these investiga- 
tions lacking the rigour, and often also the generality, desirable in view of the significance 
of the result. What is much more restrictive is that these proofs could be established only 
for energies below the three-body dissociation threshold. 

In the present paper we give, based on the three-body integral equations approach in 
momentum space, a non-perturbative derivation of the behavior of the optical potential 
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as the momentum transfer goes to zero, which is valid both below and above the three- 
body dissociation threshold, and which does not suffer from the shortcomings of the 
previous attempts. Of course, this momentum-transfer behavior determines the large- 
distance behavior in coordinate space. At the same time we deduce the general (energy- 
dependent) form of the induced dipole polarisability. For this purpose, use is made of the 
recently developed asymptotic form of the wave function for three charged particles in 
the continuum or equivalently of the asymptotic expression for the three-charged-particle 
Green function that applies to just that region of the three-body configuration space 
which is relevant in the present context. 

The paper is organized as follows. In Sec. [Ij we briefly recapitulate the formulation 
of the optical potential for elastic scattering from the m-th bound state (either ground 
or excited state) of the composite particle within the framework of the three-body the- 
ory. When writing down explicitly its leading term in the limit of vanishing momentum 
transfer, the above mentioned asymptotic expression for the three-charged-particle Green 
function is introduced. The behavior of the optical potential as the momentum transfer 



goes to zero is derived in Sec. |TTT| . There we also present the general expression of the 
induced dipole polarisability. These results are then specialized to on-shell scattering, 
both in momentum and in coordinate space. The Summary contains a discussion of our 
achievements. Several auxiliary results are collected in the Appendices. 

We mention that throughout we use natural units, i.e., h = c = 1. Furthermore, unit 
vectors will be denoted by a hat, i. e., a = a/a. 
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II. THE OPTICAL POTENTIAL 



A. Definition of the optical potential 

Consider three distinguishable particles with masses m u and charges e v , v — 1, 2, 3. We 
use Jacobi coordinates: p a (r a ) is the relative momentum (coordinate) between particles 
(5 and 7, and fi a = m^m^ / {mp + m 7 ) their reduced mass; q a (p a ) denotes the relative 
momentum (coordinate) between particle a and the center of mass of the pair (fij), the 
corresponding reduced mass being defined as M a = m a (mp + m^)j{m a + mp + m 7 ). 

The Hamiltonian of the three-body system is 

3 

H = H Q + V = H Q + Y,V V , (1) 
with H being the free three-body Hamiltonian, and 

Va = V* + V? (2) 

the full interaction between particles (5 and 7, consisting of a short-range (V^f) and a 
Coulombic part, 

= (3) 

The three-body transition operator U aa (z) which describes elastic and inelastic scat- 
tering in channel a, satisfies the equation 

U aa (z) = V a + V a G a (z)U aa (z). (4) 

Herein, 

% = J2 h a v v = v a s + v° = J2 K a v v s + £ &**y? (s) 

V V V 

is the channel interaction and 

G a (z) = (z- H a )~ l = (z-H - Va)' 1 (6) 

the channel Green function. The conventional notation 5 ua = 1 — 5 ua for the anti- 
Kronecker symbol has been used. 



Let \ipam) be the (normalized) bound state wave function (belonging to the binding 
energy E am , which we assume to be non-degenerate) of the pair (Pj) m , where the index 
m denotes the complete set of quantum numbers. The notation is to indicate that we 
allow the pair to be either in the ground or in some excited state. The free motion of 
particle a relative to the center of mass of (/3j) m is described by the plane wave |q a ). 
Then the quantity 

q Q ; E) = (c{ a \T am , am (E + iO)\q a ), (7) 

with 

/ ^am,«m(^) (Vara . ^ao I Vam) : (8) 

is on the energy shell, i.e., for 

E = E am = q 2 a /2M a + E am , q' a = q a = q ai (9) 
the physical amplitude amplitude for elastic scattering of particle a off the bound state 

Equation @ by itself does not yet lead to the desired Lippmann-Schwinger (LS)-type 
equation for the effective-two-body elastic scattering operator T am ^ am (z) since the spectral 
decomposition of G a (z) contains contributions not only from the m-th but also from all 
the other bound states, and in particular also from the continuum states, in subsystem a. 
This goal is achieved, e.g., by means of the Feshbach projection operator technique |25 
Introduce the projector onto the target state and its orthogonal complement, 



Pam | Iparn) i^am \ j Qarn 1 Pain- (-^) 



Then 



with 



G a (z) — P am G a (z) +QamG a (z) 

= G^z)+G^( Z ), (11) 



<<£|G£(*)|q*> = \^am) ^L^A i^ml (12) 

z - qi/2M a - E a 



represents a splitting of the channel Green function into a term ([12]) which is separable 
(with respect to the variables internal to subsystem a), and a remainder. Introducing the 
decomposition fllT] ) in (|) yields, e.g. through application of the AGS reduction procedure 



1~am,am{ z ) ^am,am ( Z ) ^ara,am ( Z ) „ Q /ri , , , A ^am,am(^)- (13) 

2 - Qi/2M Q + £ a 



For convenience we have introduced the relative momentum operator Q Q whose eigenvalue 
is q a . Equation ([13]) is the desired one-channel operator LS equation for the elastic 
scattering amplitude. The plane wave matrix elements of potential operator 

V^ )Qm (q' a , q Q ; Z) = (^ a \(^ am \V a + VaQam _ rr _\ QamVa | |q a ) 

Z fl a am* a^c am 

= (qLKVwlK, + v a G Q (z)v a \^ am )\q a ) (14) 

are seen to coincide with the standard definition of the optical potential. Here, the identity 

Qam\z H a QamVaQam\ Qam QamG(^z)Q a m • G (^), (1^) 

G(z) = (z — H)~ l being the resolvent of the full Hamiltonian H, has been used. Of course, 
this expression for the optical potential can also be derived directly from the Schrodinger 
equation. Its relation to the effective potentials introduced in the three-body theory has 
been described in 

The question of the solvability of an equation like (|I5D depends crucially on the singular 
behavior of the effective potential V^ am (c^ a , q a ; E) in the limit that the momentum 
transfer 

A « = q« - q« (16) 

goes to zero. In leading order the latter is defined entirely by the Coulombic part 
of the channel interaction, i.e., it does depend neither on the short-ranged part nor 
on the internal interaction V a . Of course, a given behavior of the optical potential for 
A a — > reflects itself in a corresponding asymptotic behavior in coordinate space. 
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B. Explicit expression 



In this section we establish the explicit expression of the optical potential which will 
be used in the following investigation of its analytical behavior in the limit A a — > 0. 
Let us rewrite ([H]) as 



'^am,«m( C la) Qoi] E) — V^^ am (q a , q Q ) + ^om^arrJS^ai Qai E) , (17) 



with 



K^alK, q«) = (cL\(^am\V a \lP am )\ci a ), (18) 

Vt£L{<,^E) = K\^ am \V a G Q (E + tO)V a \^ am )\q a ). (19) 



The first term (|18j) on the r.h.s. of fll7|), which is called the static potential, has for 
A a — > only the trivial Coulomb-type singular behavior 



Vam,am(qL, q«) (q' a | (Vw| Kfl^am) |<k) 

a^o 47re^ + e 7 ) + _ _ _ 

(the dots are to indicate terms ~ A" 1 and ~ lnA Q which may arise, e.g., if the target 
state is not spherically symmetric). It can be taken care of in the LS equation (O) in the 
usual manner (see, e.g., 0). 

Consequently the decisive question concerns the behavior of V^am^sCxi Qcn E) in the 
limit of zero momentum transfer. In lowest order perturbation theory, which in the present 
language is equivalent to approximating in ( |19"D the full three-body Green function G(z) by 
its lowest order approximation G a (z), and by taking into account that only the Coulombic 
part V£ of the channel potential contributes to the leading singular behavior, one has on 
the energy shell the well-known result 

V°i ( l(q' a , q«; E am ) Aa ~ ° ( q ' a \(ij am \V^G2(E am + iO)V° |Vw)|q a > Aa ~ ° A*, (21) 

or in coordinate space 

V apt(2) , \ (22) 

Vam^ctmKHa) ~ 94' V / 



for all energies. The factor a represents the polarizability of the composite particle. 
However, it is of outstanding importance to know whether this fundamental result holds 
also for the exact optical potential (pUp, i.e., even after all terms of the perturbation 
expansion of G(z) are summed up. In fact, for energies below the three-body dissociation 
threshold, i.e. for E < 0, it has been suggested in |T8|-|2"5[ that ( |2~T| ) remains true also 



for the exact expression (|T9|). It is the purpose of the present investigation to derive its 
behavior for arbitrary E, that is in particular also for E > 0. 

Let us write down the matrix element (|I9|) explicitly in the coordinate-space represen- 
tation, 

VS2L(qk q*; £) = £ £ ^ / d P ' a dv' a d Pa dr a e -*w- 

xr am K)V»(e ail p> a - \ fl r' a )G®(p' a , r' a , p a , v a , E + i0) 
xK(e Q ,p Q - A l/ r a )^ om (r Q )e- iq « ^. (23) 

For convenience, we have introduced the antisymmetric symbol €p a = —e a p, with e a p = +1 
for (a, 0) being a cyclic ordering of (1, 2, 3), and the mass ratios X u = m v j (mp + m 7 ), for 
v = /3,'y. 

To further evaluate expression (|23|) , the spectral decomposition of the kernel of G®(E+ 
iO) must be inserted and then the integrations over r a , p a , r' a and p' a must be carried 
out. This requires in principle the knowledge of all the solutions of the full three-body 
Schrodinger equation 



|E - T ra - T Pa - £ K} y (+) (r a , p a ) = 0, (24) 

whether they describe states with asymptotically two or three unbound particles. Only 
the states corresponding to the discrete spectrum of H, i.e. three-body bound states, are 
not needed since they do not contribute to a singular behavior of V°^^ in this limit (the 
reason being simply that they yield only terms which are separable in the incoming and 
outgoing momenta and thus do not depend on the momentum transfer at all). Here, T ra 
is the kinetic energy operator for the relative motion of particles (3 and 7, and Tp the 
one for the motion of particle a relative to the center of mass of the pair (J3j) m - They 
are defined as 
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.^Ei T = -^a. (25) 
2// a ' P " 2M a V ; 



The solutions of ([24]) are not known, but also not really needed. For, we are interested 
only in the behavior of V^^m f° r - > 0. Below it will be shown that in this limit the 
leading, A a -dependent term is defined by the singularities resulting from the divergence 
of the integrals in (p3|) over p a for p a — > oo, and over p' Q for — > oo. (Compare the 



discussion for an analogous two-body problem in ref. The integrals over r a and 

do not induce a singular behavior. For, the wave functions ip am (r a ) and ^£m( r a) °f the 
incoming and outgoing bound pair (Pj) m j respectively, practically confine the region of 
integration over the magnitudes of the subsystem-internal variables r a and r' Q to values 
r a,r' a < Kami where K am = ^2fi a \E am \. For this reason, when investigating the behavior 
of the leading term of in the limit A a — > it suffices to know the Green function 

G Q (p' a , r' a ; p a , r a ; E + iO) in the asymptotic region 

fi Q n Q' a , (26) 

with 

tt a : r a /p a 0, p a -> oo, (27) 

p'a -oo. (28) 

In fi„ the Hamiltonian takes the form 



H ^ H» = T ra + T Pa + Vaiva) + v c a (pa), (29) 



where 



*£(pj = lim „ {Vp(e a pp a ~ A/3r Q ) + ^ 7 (e„ 7 p Q - A 7 r a )} = — ?- (30) 

pa-*OO,r a /Pa—>0 p a 

is the Coulomb potential between the charge e a of particle a and the total charge (ep + ey) 
of the particles /3 and 7 concentrated in their center of mass. Because of this property 
v a(Pa) is termed 'center-of-mass Coulomb potential for channel a\ Note that in the 
region Q a we could neglect in fl3"0|) the short-range interactions between particles a and 
7, and a and /3, completely, and terms ~ r a / p a in the corresponding Coulombian parts. 
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The well-known class of solutions of the asymptotic Schrodinger equation 

{E-H?}9 i £Z(j> a ,T a ) = 0, (31) 
belonging to the three-body energy E = q^ 2 /2M a + E am , is given by 

*£« (A» r «) = Vw(r«)t^W (32) 

They consist of a product of the bound state wave function of the pair (Pj) m introduced 
above, which satisfies the two-body Schrodinger equation (we use the same symbols for 
operators acting in the two- and in the three-particle space) 

{E am - T ra - V a (r a )} i} am {r a ) = 0, (33) 



and the 'center-of-mass Coulomb scattering wave function', satisfying 

,0 2 



(34) 



[2M a 

The solution of (|34|) is explicitly known, 

1>$\p a ) = e*&-P.NZF{-itf a , 1; i{ql Pa - q° • pj), (35) 

with 

Here, T(x) denotes the Gamma function, and F(a,b;x) the confluent hypergeometric 
function P7| . Note that eventual bound state solutions of (53) for an attractive center- 



of-mass Coulomb potential are of no interest in the present context since they would 
correspond to a situation with all three particles being confined (see below). 

We also need the solution of the asymptotic Schrodinger equation fl3"T|) with all three 
particles asymptotically in the continuum, i.e., the one which belongs to the three-body 
energy E = g° 2 /2M a + /c° 2 /2/z a . It will be called asymptotic continuum solution of the 
Schrodinger equation in Q a . This solution has been found in p4}| , and has the form 



^S&CPa, r«) = Vffin ^(r«) e i «fi-P«. e i ^' n ^- e -^-P-) e < fa (^-^-P-). (37) 
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The Coulomb parameters are given as 



uo > 'h uo 

K f3 S 





la i 



and the momenta k° and k° are defined as the usual linear combinations of k° and q° 
\fi = e m ,n a <&/M v -\X l , with z/ = /3, 7 . (39) 



The wave function V'kojp )( r ") * s the continuum solution of the two-body-like Schrodinger 
equation 



2/j, a 



-T ra -V a (r a )^$ (p Jr a ) = 0, 



(40) 



with 



k^J = k° + ^H (41) 

Pa 

aaOJ = - E S m r&X , eavK :\ , (42) 

describing the relative motion of particles /3 and 7 with /oca/ energy E a (p a ) = 
k^ 2 (p a )/2fi a . Thus it is, in fact, a three-body wave function, the influence of the presence 
of the third particle a however being confined to a shift of the two-body relative momen- 
tum of particles (3 and 7 from its asymptotic (for p a — > 00, i.e. particle a is infinitely far 
apart) value k° to the local value k°(p a ). The latter depends explicitly on the distance 
of particle a from the center of mass of the pair (/3 7 ) m as a result of long-ranged three- 
body correlations. Nevertheless, it is to be noted that, since p a enters the Schrodinger 
equation (fHf) only parametrically via the local energy, its solutions V'ko jp )( r «) an d the 
genuine two-body bound state wave functions ?/> Q?T) ,(r Q ), which are solutions of ([33|), are 
eigenfunctions of the same Hamiltonian {T Ta + V a (r a )} to different eigenvalues, and are 
therefore orthogonal. 

Let us add two comments. 

1) It has to be pointed out that fl3~7|) is valid in all of Q a , except for the so-called 
singular directions which are defined by the conditions e aj gk° • p a = 1 and e a7 k° ■ p a = 1. 
As was shown in EM, in this whole region ^Ji£o(p a , r a ) (i) satisfies the asymptotic 
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Schrodinger equation fl3lD , or equivalently the three-body Schrodinger equation (pif) in 
Q a , up to terms of the order 0(1/ p 2 ^, viz., 

{E - H?) *$£(p a , r„) = O ( £) , (43) 

and (ii) represents the leading term in the asymptotic expansion in f2 a of the full wave 
function (p , r a ), which belongs to three asymptotically free particles in the con- 

tinuum with Jacobi momenta k° and q°, viz., 

(p«, *«) = *gg(p a , Pa) + O f-V (44) 

\P0L 



The next-to-leading term in the expansion (f44|) satisfies (^3|) in the next orders. We 
remark that the second term in (|44]) , which is of the order 0(1/ p a ), includes rescattering 
of the particles (3 and 7 from particle a Q. Its existence can easily be deduced within 
the solvable model considered in (which corresponds to = 0, m a = 00). 

2) The asymptotic Schrodinger equation (^1|) admits also another (exact) solution, 
belonging to the same energy E = g° 2 /2M a + k^ 2 /2p a as ([37]), viz. 

(P«, r a ) = ^\r«)^\p a ), (45) 
where ?/>j^(r a ) is again continuum solution of the two-particle Schrodinger equation (Hp, 



but to the energy /2p a . However, as shown in [24| this eigenfunction of the asymptotic 



Hamiltonian H® s is not the leading term in the asymptotic expansion in Q a of the solution 



of the original Schrodinger equation (24), that is, it does not satisfy a relation like 



Hence it is not to be used in the present context (as was done mistakenly in |TJ| ) . 

Consequently, the asymptotic form of the spectral representation of the three-body 
Green function G®(z), which is valid in Q a H Q' a and is therefore appropriate for the 
investigation of the singularity structure of V^^m) nas the f° rm 

GV a , r' a , p a , r a , E + zO) G^ as (p' a , r' a , p a , r a , E + zO) = 

rfqQ f ^n(r / J<t ) (P^)<o )>K (pJCn(^) 
( 27r ) 3 [E + tO- q° Q 2 /2M a - E an ] 

, / <&° ^( P ^^T( Pa ,r aj 
/ ( 27r )3 [E + iO- qf/2M a - kf /2p. W ' 
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One last point concerns the fact that in the second term on the r.h.s. of ([46]) the 
integrations over q° and k° extend also over the singular directions, i.e., directions such 
that eafik® ■ p a — 1 and e a7 k° • p a — 1 holds. But there, as mentioned above, the 
asymptotic expansion ([44]) is not valid. Thus, in principle in the vicinity of these directions 
the exact expression of the three-body wave function ^^oL (p a , r Q ) must be used, which 

1<> (V 

is however unknown. Since the latter is normalized to <5(q° — q°')<5(k° — k°'), it has 
to be integrable everywhere, including the singular directions. The consequence is that 
the contribution from the infinitely small regions containing the singular directions is 
infinitely small f28fl . Thus, when looking for the asymptotic behavior of V°^^„(q^, q^; E) 
for A Q — > 0, in the main order we can approximate the exact wave function ^„o k0 (p a , r a ) 

If) ( |; 

by ^tfiuo(Pa) r a) m the whole integration region, including the singular directions. Of 

1( v en 

course, some care must be taken before using expression ( fiT\) in the singular directions 
since, as follows from its definition (PS|), a a (p a ) diverges there. One possible remedy 
has been proposed in the Appendix of ref. |24|. It consists in writing a a (p a )/p a w 
iJ2u=t3,~/ e au^v^ p a lnF(—ii] u , 1; i(k u p a — e au k u ■ p a )). The right-hand side has the virtue 
of remaining regular when the singular directions are approached, while away from them 
it coincides for p a — > oo in leading order with the right-hand side of (0). 

Summarizing, higher order terms in the expansion ([44]) will not contribute to the 
leading term of V°^^(q' Q , q a ; E) in the limit A a — > 0, because its behavior is defined 
by the divergence of the integrals over p' a and p a in fl23|) for p' a , p a — ■> oo. And the 
contribution from the infinitely small neighbourhoods of the singular directions is infinitely 
small. Thus we can use the representation ( ^6]) everywhere. 

As discussed above, in fl a the potentials Vg and V 1 occurring in ( p3|) can in the main 
order in l/p a be approximated by their Coulombic parts Vp and . The same holds 
true in Vt' a . Thus, in the limit A Q — > one gets for the leading, A Q -dependent part V°£^m 
of the optical potential term V°^^: 

yo P t(2) , i ■ E) A< &° V opt{as) (a' a • E) 

= y^Svcfiyn / d P' a dr 'a / dp a dr a 
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(27T) 



E 



^ m [E + tO - q° a z /2M a - E an 



(2nr[E + i 0-qf/2M a -kf/2 f i c 
xV u c (e auPa - A,r a )^ m (r a )e^-^. (47) 

The radius A, which defines the lower limits in the integrals over the magnitudes of p' a 
and p a , has to be chosen so large that A ^> l/n am , with n am defined previously. This 



condition ensures that we are allowed to use the asymptotic expansion (|44] ) and to replace 
the full two-body potentials by their Coulombic parts. These approximations have already 
been introduced in (|47|). 
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III. BEHAVIOR OF THE OPTICAL POTENTIAL IN THE LIMIT OF 
VANISHING MOMENTUM TRANSFER 



In this section we will make use of techniques and results developed for an analogous 
problem in two-particle scattering in ||28|| . In Q a we can use the asymptotic expansion 

V u c (e auPa -\ u r a ) = ^ + e au \^(p a -r a ) + o(^), with v±v = P, 1 . (48) 

pa pi \pIJ 

For the following it proves useful to apply the familiar screening to the Coulomb po- 
tentials which we choose to be of exponential type in coordinate space, i.e., V£ (r a ) — > 
V£ Xr a )e~ Sra , 5 > 0. The zero-screening limit 5 — > +0 will then be taken at the end. 



Substitute (|37|) and fl48|) into (f47| ). The orthogonality of the bound state wave function 
^ami^a) to all the other bound state wave functions ip an (r a ) for n ^ m and to the scatter- 
ing wave functions ipip\p )( r a) causes the vanishing of the contribution of the first term 
of the expansion fl4"8p. Thus one obtains in leading order for A a — ► 0: 

da r r e -*(q«-<fi)-Pa . , 
^ . T^sl / dp'J dpj r 2 

<5-+0 J (27r) 3 J P ' a >A ^ a Jp a >A Ha 



V2£K(<4 q«; e) Aa ~ iim / / dp'J dp J "~ n V 



^ m [E + tO- q° a 2 /2M a - E a 

+ y rfk " 3 JT e i^ln(fcOp^-kO-p' Q ) JT e -i»jg ln(*gp„-k0.p a ) 



+ ^0 - gf/2M Q - fc° 2 /2pj J 6 p 2 



X r„ : n 2 /. a ; ,n 2 7" i ? e ^ 15 > ( 49 ) 

where we have introduced 



D nm {p a ) = e a pe a {X p e 1 - A 7 e /3 ) / dr Q ^a„( r a)(p a • r a )Vw(r Q ), (50) 
Dkg(pj m (p a ) = e af3 e a (\[3e y - \ep) J dvj^ p ^{v a )(p 0t ■ r a )^ am (r a ), (51) 

DmniPa) = D nm(Pa)i Ank° (pj (pj = D k«(p JmOJ '■ (52) 



For convenience, the plane wave has been extracted in fl49Q from the center-of-mass 
Coulomb wave function fl35|) , by writing ip^o (p a ) = e tn "'P a ^rfl\Pa)- 

in in 

Inspection of (|^) reveals that the integrals over p' a and p a will diverge for p' a , p a — > oo 
if and only if (q^ — q°) 2 — > — 5 2 and (q Q — q°) 2 — ► — 5 2 , respectively. These divergences 
will give rise to singularites in the integrand of the final q°-integration at the positions 
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qa-q°) 2 + ^ = and (q a -q°) 2 + 5 



0. 



(53) 



For 5 —>■ +0 we get from ([53]) for the positions of the singularities 



q ( 



q° = q Q - 



(54) 



Proceeding then in close analogy to |25[, it will be shown below that the coincidence of 
these singularities when integrating over q° gives rise to a sing ularity of V°^ ( ^ in the 
limit A„ — > of the form 



•\jopt (as) 

All /Lilt 



am, am \T-a> Ma; 



q a • q«- 



(55) 



That is, V°^am has a singularity in the forward direction at z — 1 in the z-plane which 
defines the behavior of the leading term in the limit z — > 1, or equivalently A Q — > 0. 

From ([54]) we can conclude that in the limit A a — > the main contribution to the 
integral over q° comes from the neighbourhood 



q° « q« ~ q Q 



(56) 



In |ZS| this has been proved in the case of two-body scattering, but it holds also for the 
three-body case considered here. For, the additional integration over k° appearing in 
the intermediate three-body continuum state contribution in ( |4^| ) does not influence this 
conclusion. Thus the leading singular behavior of V^^m ^ s defined by the divergence of 
the integrals over p a , over p' a , and over q°. 



A. Contribution from two-fragment intermediate states 

Let us first consider that part of (f49|) which is due to the contribution from all sub- 
system bound states n 7^ m in the intermediate state, 

v opt(as) . E ) A «~° lim / / do > f dn e~<-^ ^ (+) 

Km,am( b ){^^,N ~ J ^ J ^ d Pa ] ^ d Pa ^ e ^0 (pj 

x V D mn(P a )Dnm(P a ) ,(+)«, N p -fr„ g^ ^ 
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Since the behavior of V^^amib) f° r — > is determined by momenta characterized by 



, the propagator 



dJq°" 



E + iO-q a 2 /2M a -E an ' (5* 



-i 



can be taken out from under the integral over q° at the point q^ — Qa, provided it is not 
singular there. Let us write fl58|) as 

d n (q° a ) = 2M a /[q 2 an -q a 2 + iO], (59) 

with (recall (D) 

q 2 an = 2M a (E-E an ) = q 2 a + 2M a (E am -E an ), n ^ m. (60) 

If q° = q a , then for on-shell scattering, i.e. for q a = q a , d n (q a ) is always finite. However, 
for off-shell scattering (q a ^ q a ), d n (q a ) has for E — E an > a pole at q a = q an . This pole 
occurs both for E > and for E < 0, i.e. for energies above and below the dissociation 
threshold. 

Summarizing, for g° = q a ^ q an , d n (q^) is not singular and can therefore be taken out 
from under the integral over q° at the momentum g° = q a when trying to extract the 
behavior of V°^^ in the limit A a — ► in leading order |28[. Clearly this is no longer 
allowed at the point q a = q an where the propagator d n (q a ) has a pole. 

Thus, if q a ^ q an , we have in leading order 

v o P t(as) ( l . E) Ac^o v 2M a r r e -^-#» 

xD nm (K)e-^^ J ^ K?(p'J<f (P a )} ■ (61) 
We can use the completeness relation 

/ WIVJ^PJ} = ^(Pa - Pa) - E^(Pa)C(P Q ), (62) 

for the solutions of fl5I|) . Of course, if the center-of-mass Coulomb potential v^(p a ) 
is repulsive, the sum over all bound states k with wave functions ip K (p a ) is missing. 
Introducing ( |62j) into (|6T|) we obtain 
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,,opt(as) , / 77\ A ?C +0 1M a f e 1 a 'Pa 

Kmamib^^E) « d Pa |L>„ m (pJ | 2 . (63) 

n^m Ian ~ Qa •> Pc>A P a 

To arrive at this expression we have already taken into account that, owing to the presence 
of the factor p~ 4 in the integrand, the integral in (|63|) converges, so that the limit 5 — > +0 
and the integration over p a could be interchanged. 

In fact, the result ( |5"3"D is valid both for repulsive and attractive center-of-mass Coulomb 
potentials. For, when deriving it for the latter case, one only has to bear in mind that the 
bound state part in (0), when inserted in (|6T|), gives rise to terms in which the integrals 
over p a and p' a do not diverge for p a , p' a — > oo, due to the exponential decay of the bound 
state wave functions, and thus lead to nonsingular expressions. Furthermore, each of these 
terms is separable with respect to q a and q^; thus they can not contribute at all to the 
momentum-transfer dependence of V°^^m in the limit A a — > 0. 

In the case q a = q an , the propagator d n (q^) can not be taken out from under the 
integral over q° at the point g° = q a . This implies that at the discrete points q a = 
Qan, n 7^ Tn, which are accessible in off-shell scattering only, the limit A a — > in V^^mib) 
requires special considerations. This investigation is deferred to Appendix |A] where it is 
shown that at these discrete momenta V^mam{b) remains regular in the limit A Q — > 0. 

It is now an easy task to extract the behavior of fl63|) for A a —* 0. To this end we 
write (|50|) in the form 

D nm {p a ) = p a -B nm = — lim(D nm • V p )e- ip ^ , (64) 

Pa 



with 



D nm = e Q/3 e a (A/3e 7 - A 7 e/3) / dr a ^* an (v a )r a ^ am (r a ), (65) 



and use the symmetry property (p2[) . Then 

v a tS (6) (qa> q«; e) £ vpJP™ ■ v P 0(D„, m ■ v p ) j(u), (66) 

n^m Vara P ' P ^ U 

the integral J(u) being defined as 

J(u) = / dp Q — — = — dp a = — , 67 

J Pa >A p b a U J A p& 
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with u = p' — A Q — p. This integral is evaluated in Appendix |B|. Substituting its 



asymptotic form (|B.3|) for u — ► (i.e. p', p, and A a going to zero) into (|66|) leads to 



V££W<£, q«; e) a «=° c$A a + (A a ), (es) 

with 

CS } = " y E "T^T [l D -| 2 + I A « ■ D -| 2 ] • ( 69 ) 

a 

We mention that the action of V p > and V p onto the term ~ u 2 in the asymptotic form 



B.3|) for J(u) yields in the limit p, p' — > also a A Q -independent term. The latter 



was, however, omitted in flfigp since V^^L/ M was defined as the leading, A a -dependent 



am, am (6) 

contribution from the intermediate two-fragment states to the nontrivial part V^^m °f 
the optical potential in the limit A a — > 0. 

For on-shell scattering = q a = q a ) we have 

(?L - ?a)/ 2M « = 4™ - 4«, n ^ m, (70) 
and hence (|55|) becomes 

= -1 y - [id i 2 + ia d i 2 1 f7ii 

njtrn ^arn *- J a.n 

Thus, V^mam(b) as §i ven by Q68|) and ([71]) , depends only on the momentum transfer A Q , 
and not on the incoming or outgoing momentum q Q and q'^ separately, or on the energy. 
Consequently, on the energy shell we can write (pB|) as 



v opt(as) , , „ - n x v opt(as) /a \ 

"am.am (6) ^Q' ^a" "^a™ + 4U J ~ V am,am (6) 



/ dp a e^P«V^ ib) (p a ), (72) 



>p a >A 

where 

,,opt(as) /_ \ _ f l-^^(Pa)! 2 1 1 / 7 o\ 

\ n^m ^am ^an I ra 

describes the asymptotics of V^^mib) m the coordinate space. It is a local potential, 
representing the contribution to the polarisation potential from all the intermediate two- 
fragment states. 
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B. Contribution from intermediate three-particle continuum states 

We now proceed to investigate the contribution V^^mic) °f ^ ne intermediate three- 
particle continuum states to the optical potential, which is given by the integral part in 
the wavy brackets of (pE9D, in the limit A a — > 0: 



e i(q a -qO)-P Q 



xe 



-<5p Q 



Pa 



-^(q°;p a ,p a ), (74) 



with 



[E + iO-qf/2M a -kf/2 f i c 1 ' 



Consider the integration over q° in (0). Since, as discussed above, the behavior 
of V°^^ for A a — ^ is defined by momenta characterized by the restriction (|56"D , the 
function L(cQ ; p a , p' a ) can be taken out from under the integral over q° at the momentum 
q° = q^ provided it is nonsingular there. A singularity of L in the q°-plane, which may 
be of concern to us in the present context, can be generated in (|7Sp only by the pole of 
the propagator 

4o (<£) = [E + iO- q°J/2M a - k° a 2 /2fi a ] (76) 

plane. In fact, this pole induces in the function L(q°;p Q ,p' Q ) a singularity in the 
g°-plane at g° = \/2M a E (the positive square root is singled out by the infinitely small 
imaginary part +i0 in the propagator). It is an end point singularity which results from 
the coincidence of the propagator pole with the lower limit zero of the integration over 
k° 

Hence one has to distinguish two cases, 
(i) q a 7^ y/2M a E, i.e., the function L(q° ; p a , p' a ) is regular at the momentum q° = q Q 



(this is obvious since the propagator pole at k° a = J 2fi a (E — q^/2M a ) ^ does not 
coincide with the lower limit = of the integration in (|7^)). As a consequence, 
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L(q a ; p a , p' a ) can be taken out from under the integral over q a at q° = q a when looking 
for the behavior of V^mam(c) in the limit A a — > 0. This holds true a fortiori for the 
propagator d fe o(g°) in fl75p. For g° = g a , the latter can be written as (recall (|3|)) 

4° = [f a /2M a - q 2 a /2M a - kf/2^ a + E am + io] ~ X . (77) 

Note that on the energy shell, i.e., for q a = q a , expression (|77|) simplifies to 

4°(?a) = - [|#«m| + 2 /2/X a ] < 0, (78) 

that is, the propagator d fc o (g°) is nonsingular at q a = 1a = q a for all k a . 



(ii) g a = y/2M a E, which can happen only off the energy shell for E > 0. In this case 
L(q a ; p a , p'a) cannot be taken out from under the integral over q a in (|74|) at the point 

= Q« because it is singular there. A closer examination, which is deferred to Appendix 
0, reveals that in the limit A a — > 0, V^^^ is 0(A Q ) if the Coulomb potential is 
repulsive, and that it remains regular for an attractive Coulomb potential ■ In other 
words, the leading singular behavior is unaltered as compared to case (i). 

Summarizing, in order to isolate the leading singular part of V^^^ for q a ^ 
^2M a E, we can take in ([74]) the propagator d k o(q^), as well as all other ^-depending 
factors none of which becomes singular there, out from under the integral over q a at the 
point q a = q Q - This gives 

V opt(as) , , >£ ,V l im [jK. [ do' ( do e ~ l< ' P '\ ^ 



>< n \ z:_z H : ) e 



Pa 

)AcO (Pq )m{pg) f dq a e ici°-(P' a -pj 
[E + zO - ql/2M a - kf/2fia\ J (2tt) 3 

= lim [*Zf dp e '^' P \ -^ 

x Dmk° a (P a ){pa)DkO a (P a )m{pg) 

[E + iO-ql/2M a -kf/2 f i a Y 
where k° is given by the relation (139|) , but with q a instead of q°. Similarly, ff v is defined 
as in (|38f) , but with k° instead of k°. 

We remark that the fact that the leading term of (TJj) for A a — ■> is defined by that 
part of the integration region where p' Q = p a , was derived for an analogous two-body 
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problem in ||28|| . As is apparent, the additional integration over k°, which is present in 
the three-body case under consideration, does not invalidate that result. 

Next we observe that, since the behavior of V^mam{c) ^ or ^ a — ^ is defined by 
the behavior of the integrand in ( |79| ) for p a — > oo, we can approximate -Dk°(p )m(p a ) by 
-D k o m (p Q ), i.e. by the same expression ( |5"TD but with the local momentum k° (p a ) replaced 
by its asymptotic value k° which is attained in this limit. In fact, taking into account 
(|51|), (fp and (0) one sees that 



D K{pJm (Pa) = D Km {p a ) + O (±) . (80) 

An analogous expansion applies to -D m k° (p )(p«)- It is important to realize that (|80|) 
is valid for any > 0: for arbitrary positive-definite k° a always such a large A can be 
found that for all p a > A the second term in ( |8"0"D is infinitely small compared to the first 
one. This obviously does no longer hold true for = 0, as follows from the definition 



(|4T| ) of the local momentum. But the whole integral is not influenced by the behavior 



of the integrand in an infinitesimal vicinity of = if the latter does not possess a 



nonintegrable singularity there. This is verified in Appendix [A] for q a ^ \f1M a E. Thus 
we are justified in using the approximation ([8(]) in the whole region of integration over 
k° including the origin, thereby neglecting the terms of the order 0(l/p a ). 
Introducing (|80|) in (|7|) yields in the leading order 



^(cJ^.q^J ~ J Pa > A aPa pi J (2tt)3 [ E + l - ql/2M a - kf/2p a ] 1 * J 

(again, because of the convergence factor p~ A the limit 5 — > +0 could be performed before 
integrating over p a ). 

Let us write, in analogy to (|64|). 

D kgm (p a ) = — hm(D k o m • V p )e" ip -^, (82) 

Pa P^ u 

with 

D k o m = e Q/ ge Q (A^e 7 - A 7 e^) J dr a ip^* (r a )r a i m (r a ), (83) 

and similarly for -D m k° (p<J and D m ko (cf. (0)). Then we obtain from (|8T| ) 
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V 



opt (as) 
am, am (c) 



W n -^ Aa ~° linn / (Dmk ° ' V P , )( Dk > • V p) J ( U ) 

q -' q - ^ ~ p^oi (2vr)3 [£" + iO — g2/ 2 M a - fc° 2 /2/z a ] " 



The integral J(u) is given by ([H]). Making use of its asymptotic form ( |B.3| ) for w 
find 

M, q«; £) Aq = ° ^£ C) A Q + o (A a ), 



with 



C(c) = _Z1_ 



(84) 
we 

(85) 



ID 



k° I 



+ |A Q -D 



86) 



4 7 (2tt)3 [£ + i0 - ql/2M a - ' 
As in (|S8D, we omitted also in fl8"5"|) a A a -independent term. We point out that from 
the discussion in the Appendix [A] follows that the result fl85|) with ( |5ED is valid also for 
g Q = ^2M a E where, in contrast to appearance, Cj£' in fl86|) is not singular. 

For on-shell scattering (q' a = q a = q a ), taking into account ( [78]) the expression (|86|) 
simplifies to 



C (c) 



7T 







2 + 


A a • D k o m 


2" 


(2tt) 3 


| J- 1 am 







(87) 



That is, in the leading order V^^m^) depends neither on the incoming nor the outgoing 
momenta q Q and q' Q alone but only on the momentum transfer A Q . It is also independent 
of the energy. Hence we can write fl8"T|) as 



V. 



opt (as) / / , p 1 - n \ _ -iv 

am,am (c) I Via) Vlai -^am + = K 



opt (as) / * n, 
am,am(c)\ a > 



f dn e -i^ P aV opt(as) / 

J pa>A ^ ' 



where the local potential 



V 



opt (as) 



am. am c \" a 



(c)(Po 



\D 



k°mlP c 



(2vr)3 [1^1 + ^/2^] J Pi 



?9) 



describes the coordinate-space asymptotics of V^}^r c y It is the contribution to the 
polarisation potential from the intermediate three-body continuum states. 

Summing up the contribution (|68|) with fl69|) from the intermediate two-fragment states 
and the contribution fl85|) with (KH) from the intermediate three-particle continuum states, 
we thus have off the energy shell the final result 
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V££2(<£, q«; *0 Aa = Ci A « + ( A «X ( 9 °) 

with 

/-? _ ^7(6) I We) 




| 2 + |A Q -D nm | 2 r dk° |D k o m | 2 + |A a -D;- 12 



+ 



'k°m 



(91) 



4 [ n ^[£ + *0-g 2 /2M Q -£ an ] J {27if[E + tQ-ql/2M a -kf/2 l i a 
Consequently, for off-shell scattering we have found for the asymptotic behavior of the 
total optical potential (|17D , i.e. including the contribution from the static potential (for 
a spherically symmetric target state m), in the limit A a — ► the result 

VSU(<& E ) A =° 47re " ( A2 +e7) + ^ A « + ( A «)> ( 92 ) 

a 

with the energy- and momentum- dependent factor G\ = Ci(q a ,E) given by (pl~|). It 
should be noted that the strength factor C\ of the leading term decreases with increasing 
energy. 

Equation (j9~2"D with ( pT|) constitutes our first main result. It implies the exact compen- 
sation in the (non-static part of the) optical potential of the singular terms proportional 
to A^ 1 and lnA a , in contrast to their appearence (as discussed in |20| for negative en- 



ergies) in the effective potentials which occur in the effective-two-body formulation of 
the three-charged particle theory. We emphasize that this compensation was proved here 
both on and off the energy shell, and for arbitrary total three-body energy E < and 
E > 0. It generalizes the cancellation proofs, given in |IS| - |2"3|| for energies below the 
dissociation threshold, to positive energies. In addition, at negative three-particle energy 
our derivation avoids the various inconsistencies existing there. 

If we consider on-shell scattering then, as follows from (0), ([71]), ( |35D and (|87|), V°^^ 



depends on the momenta only in the form of the momentum transfer A a , and no longer 
on the energy, i.e., 

Consequently, the inverse Fourier transform of V°^^(A a ) is a local, energy-independent 
potential describing the asymptotic behavior of the (non-static part of the) optical po- 
tential in coordinate space. From ( |73|) and ( |89| ) we read off 
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^(pJ^-A + ofi), (94) 

Z Pa \PaJ 

with 

Q = 2E l^ m (pJI 2 +2 /j% JAc° m (pJI 2 _ (g5) 

l^aml - | -Ban | J i^V \E am \ + k J /2fi a 

As can be seen, the strength factor governing the coordinate-space behavior is just the 
static dipole polarisability as known from the perturbative approaches. That is, no re- 
mormalization of a arises from summing up all the higher order terms in the perturbation 
expansion of the full three-particle Green function G(z). Furthermore, all dependence on 
the incoming energy has disappeared (in contrast to the approximate result of ref. [I7| , 
but in agreement with the result obtained in ]nj within the adiabatic approach). 



For the coordinate-space behavior of the total optical potential, i.e. the inverse Fourier 
transform of fll7f) , we derive on the energy shell for spherically symmetric targets 



V opt r x pc^oo e«(e/3 + e 7 ) _ _a_ f 1_ i , , 

y am,am\r'aJ rj 4 1 1 4 I \ UKJ J 

Pa Z Pa \ Pa , 



Equation ( p6|) with (p5|) is our second main result. It states that in the asymptotic 
expansion of the energy-shell restriction of the optical potential for large intercluster 
separation there occurs as the first nonvanishing term after the center-of-mass Coulomb 
interaction, which arises from the multipole expansion of the folded Coulomb channel 
interaction, the local (induced) polarisation potential —a/2p^. This result extends the 
one derived in |T8|-p^1 for E < 0, to arbitrary energies. 

Simultaneously we derived the exact expression for the static dipole polarisability a. 
To our best knowledge this is the first derivation of a within the framework of a genuinly 
three-body, non-perturbative, non-adiabatic approach. In this context we point out that 
the same expression for a was derived in jn| for E < 0. However, in that derivation a 
spectral decomposition of the three-body Green function in Q a H VL' a was used, which is 
inadequate since it does not take into account the long-ranged three-body correlations 



described in [24]. The latter necessarily lead to the appearence of the wave functions 



■0kt(p )( r a) instead of ^(r Q ). Though, this inconsistency in [18| did not influence the 
final result. Nevertheless, we feel that it is important to present a consistent three-body 
derivation of the dipole polarisability. 
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IV. SUMMARY 



We have investigated the singularity structure in momentum space of the optical po- 
tential responsible for the elastic scattering of a charged particle a off a bound state of 
two charged (or one charged and one neutral) particles, on the basis of the rigourous 
three-body scattering theory. Without having recourse to perturbative or other approx- 
imate methods we have shown that (for spherically symmetric bound states) the op- 
tical potential behaves, in the limit that the momentum transfer A a goes to zero, as 
C A~ 2 + CiA a + o (A Q ), for both negative and positive energies. That is, terms ~ A~ x or 
~ In A a cancel exactly. The factors appearing in this expansion have the familiar physical 
interpretation: C is the product of the total charges of the colliding fragments, and C\ 
is the general, energy- and momentum-dependent expression for the polarisability of the 
composite particle. 

For on-shell scattering we find that this momentum-transfer behavior entails in co- 
ordinate space a local tail of the optical potential of the form C /p a — a/2p A a + o(p~ 4 ), 
where p a denotes the distance of the elementary from the center of mass of the composite 
particle. The polarisability a derived here coincides with the expression as extracted in 
perturbative or similar approximate approaches. These results imply that 

(i) the exact cancellation of terms ~ p~ 2 and ~ p~ 3 takes place for energies below and 
above the dissociation threshold, 

(ii) no renormalization of the polarisability a due to summing up the infinite perturbation 
series occurs, and 

(iii) the strength factor a which governs the long-distance behavior in coordinate space is 
independent of the incident energy (in contrast to its off-shell analog C\ which determines 
the strength of the corresponding singularity in momentum space). 

t Permanent address: Institute for Nuclear Physics, Tashkent, Usbekistan 
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APPENDIX A: 



In this Appendix we first investigate the behavior in the limit A Q = |q^ — q a | — > of 
the contribution of the n-th intermediate bound state in expression (|57D , to be denoted 



by ^amam(n)> a ^ the discrete point q a = q an . To simplify the considerations we replace 
the Coulomb scattering wave functions ^io^*(p Q ) and ip^o\p' a ) by their leading terms for 

io Ma 

large distances, namely by the planes waves e _iq « Pa and e l<i °'P^, respectively. We also 
omit the functions D nm (p a ) and D mn (p' a ). By doing so the question of the integrability or 
non-integrability of the singularities under considerations is not influenced. We therefore 
consider 

V££L>(<4 q«; E) 2M a Km ^ _ dp' a / dpf 



6^+0 J [Any Jp' a >A Jp a >A p, 
1 c - e *(«fi-«ia)-Pa 

x — r^e-^ 



A a ^0 

~ iim 



[q 2 an + iO-q° a z ] Pi 
dqj 1 1 

02l 



<^+° J W v/(q° - q a ) 2 + <5 2 [<&» + 

1 



(A.l) 



(<fi -q a ) 2 + 5 2 

The integrals over p a and p' a from to A added in the second line are certainly less 
singular than the original expression; hence they do not alter our conclusion. We use the 
integral representation 

Tchf^W-clC^TvW-C (A ' 2) 

where § means the integral along a closed contour around the point (3 = in the 
integral from to oo of the second representation we must have arg (3 = 0. This leads to 

V££2 („)(<& q- E) Aq ~° Km jT d/3-j= f da^=J((3, a; q' a , q Q , q an ), (A.3) 



with 



ion J 



J(P,a;c^ a ,q a ,q 

f dq° a 1 1 1 

7 (2vr)3 \p + (qo - q ^)2 + p] [q 2 n + l0 _ g o2] [a + (q o _ ^2 + ^2] • 

This integral is explicitly known | ]29| | , 

29 



(A.4) 



o7r a/c 2 — ao 



c + a/c 2 — ab 
c — \Jc 2 — ab 



(A.5) 



with 



ab = [A 2 a + (v^+ V^) 2 M + (v^ - ^n) 2 ][gl 2 + (V^ - V) 2 ], (A.6) 
c = -zg an [A* + (v^ + Jftf] + y^fg 2 + a _ g 2 j + + ^ _ ? 2j (A>7) 

Because of the presence of the parameters a and ft, the regularisation parameter 5 is 
no longer needed. We therefore have already put 5 equal to zero in ( |A.5|) - (|A.7|) . For 
Qa = q' a = Qan, flOD and (TAT7D simplify to 



ab = ^^ft[A 2 a + (V5 + ,/5) 2 ] [v^ - 2ig an ] - 2ig w ], (A.8) 
c=-ig an [A 2 + ( v / ^+ x /(3) 2 ] + yfi/3[y/Ei + yf(3]. (A.9) 



Substitution of a = A^w and ft = A 2 a t in ( |A.8|) and ( |A.9|) shows that in the limit 
A a ^0 



a6~A*, c~A 2 , i.e., Vc 2 - a& ~ A 2 . (A.10) 



Introducing ( A.5 ) into (|A.3| ), changing to the new variables v and t, and taking into 



account ( |A.10| ) , it is immediately seen that the main term of ^ , n s is independent of 
A a in the limit A Q — > at the discrete point g a = g an , i.e., it is regular there. Hence the 
same holds true also for the sum over all intermediate bound state contributions V^^m. 
The other problem concerns the behavior of the intermediate three-body continuum 



part @ if 

Qa = Qai, where q al = ^2M a E. In that case we are not a priori allowed to 
take out the propagator (ifco(g°) from under the integral over q° at the point g° = q a . 
To investigate this case consider the definition ( [75] ) of L(cQ; p a , p' a ) . Since in the limit 
A Q — ► the leading, A a -dependent term is defined by p a « p' a , the exponential Coulomb 
distorsion factors cancel. That is, it suffices to consider 

Hal ; Pa , P ' a ) - / D ^^^P^ (A.11) 

If we perform the k°-integration, the pole of dko(q^) can give rise to a singularity 
of £(q°; p a , p' a ) in the g°-plane at g° = q a ±. It is an end point singularity which arises 
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from the coincidence of the propagator pole at k° a = \J fi a (Q a i ~ Qa 2 )/M a with the lower 
limit k\ = of integration over Since the leading, A Q -dependent term of V°^ ^^ 
in the limit A Q — > is generated by the coincidence of the singularities of the integrand 
in the integral ([74]) over q° at the point q° = q a = q a , the appearence of an additional 
singularity, namely in the the function L, will influence the leading term of V°^^*( c \. 

To investigate the analytic behavior of L(cQ; p a , p' a ) for g° — g al — * 0, consider the 
integrand in ( |A.11| ). The denominator has a zero that can cause the singularity under 
consideration. The behavior for k^, — > of the numerator (cf. (|51|)) is defined by the wave 
function iptf)\p )( r a) which we write as 

1 



^Wj( r «)= ex P 



a a{P a ) _ 

V 1,0 



Pa 



- ^(r a ) + ^'^ ^(rj +o(\). (A.12) 

The threshold behavior for k^ — * of the scattering wave function ^>£p (r a ) in the potential 
Va — Va + K* * s entirely given by the corresponding behavior of the Coulomb scattering 
wave function V'rwo ( r a)- This follows from the wellknown fact that the threshold behavior 
of the radial wave functions for two charged particles is the same for all the partial waves 



and is defined by the factor |30| 

2^ 



, „ , . (A.13) 

Therefore, the behavior of the full wave function ip^t\r a ) for k° — > is governed by the 
same factor, and the one of the numerator in ( |A.11| ) by | N a \ 2 . 
Hence we have to distinguish two cases. 

1. The Coulomb potential V£ is repulsive, i.e., rj^, > 0. Then for k° a — > one has 
| N a | 2 « 27T?7°e _27rT '" — > 0. The exponential smallness of the numerator for — > 
compensates any zero of the denominator at the point = 0. In other words, the pole of 



the propagator in (|A.11|) at k^ = y ^ a (q a i ~ Q a )/M a does not generate a singularity of 
Pa-, Pa) at tlie P oint Qa = Thus, when extracting the leading term of VZ!i,am(c) 
for A Q — > 0, the function L(q° ; p a , p' a ) can be taken out from under the integral over q° 
at the point q° = q Q , as in the case q a ^ q a i- The consequence is that the behavior of 
the optical potential is the same for all q a including q a = q a i, namely ~ 0(A a ). 
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2. The Coulomb potential V£ is attractive, i.e., rf^ < 0. In this case, | N a | 2 ^27r?7° — > 
oo for — > 0, and hence all terms in the expansion ( |A.12j ) become singular. This 
difficulty can be overcome as follows. According 153] the wave function \I/^tw) S (p , r a ), 
as given by (|37]), is asymptotic solution of the three-particle Schrodinger equation in 
the region Q a , i.e., it satisfies this equation up to terms 0(1/ p^) (cf. ((43[)). From the 
derivation given there it is, however, easily seen that for = 0, taking into account the 
local momentum k° (p ) in the asymptotic solution ^jij>ijf (p„, r a) instead of its asymptotic 

H.CH (V 

value k° = 0, gives corrections ~ 0(1/ p^) to the Schrodinger equation in Q a . Therefore, 
in a small vicinity of k® = the asymptotic solution ^ioi^ffPa, r a ) may be replaced by the 

1(1 Ot 

function ^(roje^^ae^^^-^'^e^ 1 "^- 6 ^-^, i.e., by (|7]) but with k° 
substituting k° (p a ) in ^ko jp )( r a)- Hence for arbitrary we may rewrite the asymptotic 
solution ^^)^o(pa) r ») in the form ( |3"7| ) with 

k°(pJ = k°+ X (0^, (A.14) 

Pa 

where is a characteristic function, which equals one everywhere except for a small 

neighbourhood of = 0, and goes smoothly to zero when /c° approaches zero. Thus, 
taking into account the factor fc° 2 from the phase volume, the behavior of the numerator 
in (|A .11 ) in the limit k° a — > is governed by k°^\ N a | 2 ~ k° a . 

Let us write £(q°; p a , p' a ) as the sum L = L e + L e , where L e denotes the integral 
over k° over the interior of a small sphere with radius e. L £ contains the remaining 
integral. Since the latter does not contain the origin /c° = 0, L £ (q^; p a , p' a ) is regular at 
Q a = lai = Qa and can therefore be taken out from under the integral over q° in ([74]) at 
the point q° = q Q . The consequence is that at the momentum q a = g al , its contribution 
to the optical potential V£ is ~ 0(A a ), i.e. of the same order as for q a ^ q al . 

However, L e (q^; p a , p' a ) is singular at g° = q a \. When performing the k°-integration, 
the pole of <4o(g°) gives rise to a behavior L £ (q°;p Q ,p^,) ~ ln(g° 2 + iO - q^) in 
the g°-plane. In order to arrive at this result we singled out from the product 
-DmkO(p^)(p' Q )-Dko (p a )m(p a ) the factor | N a | 2 which is singular at k° a = 0, and ig- 
nored the rest since it remains finite. Denote by V^^/ e \(q^, q a ; E) that part of 
K P Zmc)(<i^E) which contains L £ (c£; p a , p' a ) instead of L(c£; p a , p'J (cf. (0)). 
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Substituting ln(g° + iO - q^) for £ e (q°; p a , p' a ), and integrating over p a and p' a , gives 
v o P t{as) , , . F] A a ~,o . /■ In (g| + zO - g° 2 ) 



When writing down this expression we already took into account that the behavior of 
Kmam^K'l"'^) f° r — ^ is defined by the region p a mp' a where the product of 
the Coulomb distortion factors equals one. The integration region over q° in ( |A.15|) 
supposedly consists only of the neighbourhood of q° ~ q« ~ q„- Comparing with (|A.1|) 
(putting there q an equal to q a \) we conclude that expression (|A.15| ) is less singular at 
A a — > than ( A.lj) . Hence, in the case of an attractive Coulomb potential V^, for 
q a = q Ql the contribution y££J (<0 (<£, q Q ; E), and thus also V^^^E), is 
finite in the limit A„ — > 0. 
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APPENDIX B: 



In this Appendix we evaluate the integral ( |5"7| ) and extract its asymptotic behavior 
for u — > 0. Substituting the new variable v = up a and making use of eqs. 3.761(2) and 
8.352(3) of |JT] we obtain 



J(u) = 27rm 3 [r(-4, iuA) - T(-4, -iuA)] 



ITXU" 

~12~ 



[T(0, iuA) - r(0, -iuA)] - £ + (-1)™^) , (B.i; 



m=0 



where T(/i, x) is the incomplete Gamma function [31]. With the help of eqs. 3.761(2), 
8.230(1) and 8.232(1) of || one finds 

(B.2) 



00 ( — I^+VtM") 2 " 1-1 

r(0, iwA) - T(0, -iuA) = 2% s\(uA) = -in + ^ " 



1 a (2m- l)(2m- 1)!' 

where si(x) is the sine-integral. Substitution of (P.2j ) into (|B.1|) leads to the following 
asymptotic behavior of J(u) for u — > 0: 



..„_>() 47T 27T 2 7T 3 . 4 . 

■/(w) = — T7 H it + (J\u ). 

v ; 3A 3 3A 12 v 1 



(B.3) 
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